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This work presents the most detailed and extensive theoretical study to date of the structural 
configurations of Ruddlesden-Popper (RP) phases in antiferrodistortive (AFD) perovskites and 
formulates a program of study which can be pursued for RP phases of any AFD perovskite system. 
We systematically investigate the effects of oxygen octahedral rotations on the energies of RP phases 
in AFD perovskites (A„+iB n 03 n +i) for n = 1 ... 30, providing asymptotic results for n — > oo that 
give both the form of the interaction between stacking faults and the behavior of such stacking 
faults in isolation. We find an inverse-distance interaction between faults with a strength which 
varies by as much as a factor of two depending on the configuration of the octahedra. We find that 
the strength of this effect can be sufficient to (a) stabilize or destabilize the RP phase with respect to 
dissociation into the bulk perovskite and the bulk A-oxide and (b) affect the energy scales of the RP 
phase sufficiently to constrain the rotational states of the octahedra neighboring the stacking faults, 
even at temperatures where the octahedra in the bulk regions librate freely. Finally, we present 
evidence that the importance of the octahedral rotations can be understood in terms of changes in 
the distances between oxygen ions on opposing sides of the RP stacking faults. 

PACS numbers: 61.50.Ah,61.50.Lt,61.50.Nw,61.66.Fn,61.72.Nn,68.35.Ct,68.35.Dv,68.65.Cd,81.05.Je 



I. INTRODUCTION 

Perovskites possess a vast range of scientifically inter- 
esting and technologically important properties. These 
materials arc highly valued for their dielectricityP^ 
ferroelectricity semiconductivityj^"^ superconductiv- 
ity EMO] catalytic activityp^^ and colossal magnetore- 
sistanceP^ Such physical properties allow for the use 
of perovskites in diverse technological applications, in- 
cluding tunneling semiconductor valves and magnetic 
tunnel junctions in spintronicsp^^ dielec tric insulators 
in dynamic random access memory, 22 23 thin films in 
graded ferroelectric devices (GFDs)j- 24 * 25 ! and alterna- 
tive gates in meta l-oxide-semiconductor field transistors 
(MOSFETs)P23 

Cation stoichiometry has a large impa ct on the physics 
of these materials, but is difficultpSSD to control dur- 
ing film growth. Cation non-stoichiometries in the per- 
ovskites can potentially take the form of point or pla- 
nar defects, with perovskites experimentally found to re- 
sist point defect formation upon sufficient cation non- 
stoichiometryPSEI! i n A 2+ B 4+ Og _ materials (such as 
SrTi0 3 ) with an excess of species A (or, more pre- 
cisely, additional AO to maintain charge neutrality), 
the resulting planar defects produce a series of homol- 
ogous compounds of the form A„ +1 B„0 3n+ i. (Note that 
bulk ABO3 corresponds to one end member of this se- 
ries, n = 00, and bulk AO corresponds to the other 
end member, n = 0.) The A n+ iB n 03„+i com poun ds, 
known as the Ruddlesden-Popper (RP) phasesp322 re- 
flect the modification of the cation stoichiometry by the 
addition of an extra unit of AO per n units of bulk 
ABO3 material. Structurally, they take the form of bulk 
ABO3, separated every n layers by the insertion of ex- 



cess (001) planes of AO to create stacking faults of the 
form . . ./AO/BO2/AO/AO/BO2/AO/. . . in the normal 
. . ./AO/BO2/AO/BO2/AO/. . . stacking sequence of the 
material. Across each such stacking fault, the bulk per- 
ovskite slabs are alternately displaced by in-plane vectors 
of the form ^[±1, ±1, ]. Figure IT] depicts four members 
of the homologous series A n+ iB,703n+i, with n = 1, 2, 
3, and 00. 

RP phases, specifically, manifest additional properties 
beyond th ose o f bulk perovskitesjthese include high di- 
electricityj 35 * 36 ! thermoelectricit y photocatalysisp^ un- 
conventional superconductivity]^ quantum criticalityj^ 
metamagnetismpf ferromagnetismp2 colossal magneto- 
resistancej^ and x-ray optic adaptability. 44 Despite the 
interest in such compounds, experimentalists have faced 
difficulties growing RP phases for "intermediate" val- 
ues of n (n > 4) ! 45 " 48 l For example, strontium titanate 
(SrTi0 3 ) compounds with such stoichiometrics tend not 
to form growths of a single uniform RP phase, but rather 
materials with multiple regions, each with RP phases of 
different values of n, so-called "inter growths ,"i 49 " 5 1 1 In- 
deed, conventional ceramic powder sintering has only 
produced single-phase R P sam ples for n < 3: SraTiOi, 
Sr 3 Ti 2 7 , and Si^TiaOio- 13 ^ More recently, successful 
creation of intermediate members of the Sr n+ iTi n 03„+i 
series has been accomplished with advanced techniques 
for epitaxial growth of thin-films under thermodynamical 
non-equilibrium conditions, such as sputtering,^ pulsed 
laser de posit ion (PLD)J 53 " 55 land molecular beam epitaxy 
(MBE)pSEMSS w ith characterization by high-resolution 
transmission electron microscopy (HRTEMp-3 anc j x-ray 
diffraction (XRD)|S21 Such recent observations renew the 
urgency for a careful theoretical study of RP phases with 
intermediate and large values of n. 
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FIG. 1: Four members of series of homologous compounds of form A n +iB n 03„+i, showing n layers of bulk perovskite between 
stacking faults: (a) A2BO4 with n — 1, (b) A3B2O7 with n — 2, (c) A4B3O10 with n = 3, and (d) bulk ABO3 with n = 00. 



Previous theoretical work on RP phases of perovskites 
includes a limited number of both ab initio electronic 
structure calculations and empirical potential studies. 
A pioneering empirical-potential stud}®' in the 1980's 
considered all values of n < 12 for strontium titanate 
but did not report sufficient precision to resolve differ- 
ences among any of the phases above n = 1. A decade 
later, a subsequent empirical-potential studjElwas able 
to resolve differences among phases but only considered 
phases for n < 4. Total energy electronic structure stud- 
ies have also been carried out, primarily for strontium 
titanate^HlU ( n < 3) but also for three additional per- 
ovskite transition-metal oxides^ (n = 2 only). None of 
these studies have explored the significant interactions 
which we find to exist between the rotational states of 
the oxygen octahedra on opposite sides of the AO/ AO 
stacking faults. 

Many perovskites that form RP phases exhibit anti- 
ferrodistortive (AFD) behavior in which the oxygen oc- 
tahedra tend to rotate in an alternating spatial pattern 
from their ideal orientations with a relatively low energy 
scale, so that numerous structural phases exist for the 
corresponding bulk materials. Unless a full quantum sta- 
tistical treatment of the RP phases is considered in the 
theoretical calculations, the bulk regions in these phases 
will tend to relax to the zero-temperature ground state 
within the model used to describe the material. One 
must therefore carefully and explicitly consider two in- 



terrelated effects: (a) the interaction between different 
possible antiferrodistortive reconstructions on either side 
of the stacking faults and (b) the possibility that these re- 
constructions will assume different orientations and thus 
prefer lattice structures in potential conflict across the 
stacking fault. If the different possible reconstructions of 
octahedral rotations are not considered, the true ground 
state of the RP phases may be missed. If the lattice 
vectors are not allowed to relax fully, the extraction of 
the formation energy of the RP phases will suffer errors 
which scale with the size of the bulk regions (i.e., lin- 



early with n), as we demonstrate below in Section VA 



The aforementioned theoretical works thus have limitec 
applicability, even for the small values of n which they 
consider. 

Indeed, experiments reveal that these effects are known 
to be present in RP phases of perovskites; for example, 
both Sr n+1 Ti n 03„ + i and Ca„ + iTi„03„ + i exhibit rota- 
tions of their oxygen octahedra, for n — 2 and n — 3 RP 
phases.®! Previous work has considered general group- 
theoretical studies of the space groups of po ssible oc- 
tahedral rotation patterns in RP phases E2IZU but with 
restrictions either to the n = 1 RP phase or to a lim- 
ited number of possible relative rotational orientations 
on opposite sides of the stacking faults. 

In this work, we provide a careful, comprehensive 
atomistic study of the Ruddlesden-Popper phases of a 
physical model of an antiferrodistortive perovskite, con- 
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sidering a wide range of n (specifically, n = 1...30), 
properly accounting for the full lattice relaxation of the 
RP phases, and exploring all combinatorial possibilities 
of different orientations of oxygen octahedra on opposite 
sides of the AO/ AO stacking faults. In Section [TT] we 
introduce the class of materials which we consider, those 
antiferrodistortive perovskites which form intrinsic RP 
phases and belong to the Glazer system of the greatest 
symmetry commonly found in nature: a~a~a~y^ Sec- 



tion [HI] then presents an exhaustive catalogue of the dif- 
ferent possible structures for RP phases in perovskites 
within this system. In generating this catalogue, we in- 
troduce a convenient symmetry algebra which allows one 
to quickly enumerate the structures, and we use it to find 
a total of five symmetry-distinct possibilities (consistent 
with distorted phases Nos. 18-22 which Hatch et alP-3 
found for the n = 1 RP phase after extensive searching 
with a computer program). Section IV then introduces 



the shell-potential model and the numerical techniques 
which we use. 



Section |V A| explores the aforementioned five distinct 
configurations for all RP phases for n — 1 ... 30. We 
find that the energies of RP phases are indeed quite sen- 
sitive to octahedral rotations, sufficiently sensitive that 
unfavorable configurations become unstable relative to 
phase separation into bulk perovskite and bulk A-oxide. 
In fact, the effect is sufficiently strong to suggest some in- 
triguing possibilities. For low densities of stacking faults 
(high n), the rotational states of the octahedra neighbor- 
ing the faults might be constrained, even at temperatures 
where the octahedra in the bulk librate freely, so that dif- 
ferent degrees of order are observed in the bulk and at 
the interfaces. For high densities of stacking faults (low 
n), this effect may increase the transition temperatures 
associated with the octahedral rotations. 

Section IV Bl considers interactions between the stack- 
ing faults present in the RP phases. For each config- 
uration, we examine the energy of the A n+ iB n 3n +i 
RP phases as a function of n, which directly measures 
the separation between stacking faults. We demonstrate 
that the interaction is quite sensitive to the octahedral 
rotations, differing in strength by as much as a factor 
of two depending on the rotations. Next, Section VC 



examines the issue of the asymptotic form of this inter- 
action for tl — ► oo . We find that the interaction between 
stacking faults varies as the inverse of the distance be- 
tween them, and we extract both the binding energies 
of stacking faults and the formation energies of isolated 
faults. Again, we find that the interaction energy be- 
tween faults is highly sensitive to the different possible 
rotational states of the oxygen octahedra and may even 
lead to ordering at the stacking faults at temperatures 
where the bulk regions have lost their oricntational or- 
der. This section then concludes with a proposal for a 
simple physical mechanism to explain the strong depen- 
dence of the interfacial energy on the rotational state of 
the octahedra: some configurations result in movement 
of like-charged neighboring oxygen ions directly toward 



each other and thus are high in energy, whereas others 
result in movement of oxygen ions past each other and 
thus are low in energy. 

Finally, in Section |VD| exerting care in tracking the 
chemical potentials of the various reference systems, we 
verify that the formation of isolated planar stacking 
faults is favored over the formation of point defects in 
our model material, regardless of the rotational state of 
the oxygen octahedra. 



II. BACKGROUND 

In this work, we embark on a study of the generic ef- 
fects of octahedral rotations on RP phases. For this first 
such study, we shall focus on the most highly symmetric 
of the twenty-three possible Glazer systems PS There are 
two such systems, denoted in Glazer notation as a + a + a + 
and a~a~a~ , only the latter of which is commonly found 
in nature. Exam ples o f a~a~a~ perovskites include 
LaAlOag NdA10 3 pEl CeA10 3 P BiFe0 3 jEl LiNb0 3 JEI 
LiTa0 3 P PbZr . 9 Ti .iO 3 (PZT),ES and many others 
catalogued by Glazer^! and Megaw and DarlingtonPZl 

The Glazer notation refers to the relative state of ro- 
tation of neighboring oxygen octahedra in antiferrodis- 
tortive reconstructions of the perovskite structure. In 
actual fact, the motion of the octahedra within such re- 
constructions is only approximately a rotation as the oxy- 
gen atoms are confined to the faces of each cube; regard- 
less, we shall refer to their motion as rotational through- 
out. Glazer exhaustively catalogued all twenty-three pos- 
sible patterns ofthese octahedral rotations ( "tilts" in his 
terminology) with each category assigned an appro- 
priate nomenclature to denote the axis of rotation and 
the relative sign of successive rotations along that axis. 
In brief, in the a~~a~a~ system, all oxygen octahedra 
rotate either clockwise or counterclockwise about a fixed 
trigonal axis in an alternating checkerboard pattern in all 
three dimensions, resulting in a cell-doubling reconstruc- 
tion with a 2 x 2x2 supercell. (Greater details appear 
in Section III ) 

The antiferrodistortive phase transition, associated 
with these rotations, occurs due to the softening of 



the optical phonons at the R corner ([111] zone 
boundary) of the Brillouin zone, a s fam ously studied 
experimentalises] anc j theoretically^ 81 * 82 ! and reviewed 
extensively.^ The transition is therefore also sometimes 
known as a Zone Boundary Transition (ZBT) and was 
recognized even earlier as a means by which crystals could 
double the size of their primitive cells The rotation of 
the oxygen octahedra accompany the softening of these 
phonons, with their rotation angles serving as the order 
parameters.^ The degeneracy of this phonon mode in the 
higher-temperature cubic phase enables the rotation of 
octahedra about different axes — (100), (110), and (111) 
— to break the crystal's symmetry below its critical tem- 
perature. While the a~a~a~ system results from rota- 
tions about the (111) axis, all three systems emerge from 
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the same underlying instability J 83 * 86 ! with anharmonic in- 
teractions determining the resulting low-temperature lat- 
tice symmetry.^ As the underlying physics is so similar, 
we thus expect that perovskites in Glazer systems other 
than a~a~a~ will exhibit similar generic behaviors to the 
results expounded below. 

In seeking a model potential for our study, we searched 
the literature for shell-potential parameters among the 
a~a~a~ perovskites. Unfortunately, we discovered that 
the ground states of the available models generally do not 
correctly match their corresponding Glazer systems.' 8 - 8 -' 
We were only able to identify a single material, lan- 
thanum aluminate (La 3+ Al 3+ 03~), with a ground state 
in the correct Glazer system Unfortunately, lanthanum 
aluminate does not form intrinsic RP phases, since 
its composition as A„ + iB„03 n +i would violate basic 
charge balance; intrinsic RP formation requires per- 
ovskites with an A 2+ B 4+ 2_ chemical formula for the 
additional A 2+ 2 ~ layer to be neutral. Lanthanum alu- 
minate, however, can form extrinsic RP phases by in- 
corporation of additional neutral layers of another per- 
ovskite, strontium oxide: SrO • La„Al„03„. 

For simplicity of this initial theoretical study, we focus 
on perovskites that form intrinsic RP phases, examples of 
which from the a~a~a~ system do indeed exist in nature 
(e.g., BaThO^. However, we are not aware of shell 
potentials for any of these materials. On the other hand, 
we discovered that the shell-potential model commonly 
used for strontium titanate (Sr 2+ Ti 4+ Og~), which forms 
intrinsic RP phases, does possess a ground state of the 
a~a~a~ type. 

Indeed, on a microscopic level, strontium titanate is 
very similar to l antha num aluminate. In fact, early 
x-ray experiment d 89 * 90 -' erroneously predicted a low-tem- 
perature (T < 35 K) phase transition in strontium ti- 
tanate to a rhombohcdral lattice (generally associated 
with a~a~a~ microscopic ordering), only t o be corrected 
by subsequent spectroscopic studies .EMU Strontium ti- 
tanate assumes its actual tetragonal ground-state struc- 
ture through a famous cell-doubling antiferrodistortive 
phase transition near 105 K, 92 whose "physical origin is 
the same'^as that of lanthanum aluminate. Ind eed, th is 
transition is "strikingly analogous in all respects"' 93 ' 94 ' to 
that of lanthanum aluminate, although the energy scales 
are quite different. (The structural phase transition tem- 
perature for lanthanum aluminate is 800 K. ' 85 * 95 ') Further 
connections between these two perovskites are enabled 
by the classical concept of the Goldschmidt tolerance 
factor pSlf = i rA±m a normalized ratio of the radii of 

' V2 r B +ro ' 

ABO3 ions historically used to categorize perovskites and 
predict their ground states.^ In fact, lanthanum alumi- 
nate has a tolerance factor within 1% of that of strontium 
titanate (data compiled by Shannon^). 

In sum, the parameters for the shell potential com- 
monly used for strontium titanate describe a perovskite 
which intrinsically forms RP phases within the desired 
a~a~a~ Glazer system. It is not surprising, then, to 
discover empirical evidence for perovskites, such as lan- 



thanum aluminate, which closely resemble strontium ti- 
tanate but with a~a~a~ octahedral rotations. The stan- 
dard shell potential for strontium titanate can thus rea- 
sonably be utilized as a parameterization for generic 
a~a~a~ perovskites which form RP phases. Therefore, 
in interpreting our results, one should be mindful that 
details, such as the values which we find for the vari- 
ous quantities, may not apply to any specific perovskite. 
However, the general phenomena that we uncover, such 
as the enumeration of possible octahedral configurations 
and the general form and scale of the various interactions, 
may be taken as representative of the class of a~a~a~ 
perovskites in RP phases. Moreover, since the underlying 
physical mechanism is the same for other Glazer systems, 
many of the phenomena which we uncover should be con- 
sidered for the RP phases of all AFD perovskites. 



III. CLASSIFICATION OF OCTAHEDRAL 
ROTATIONS IN RP PHASES 

We return now to a careful examination of the specifics 
of the a~a~a~ Glazer reconstruction, as visualized in 
Figure [2] The oxygen octahedra each rotate around 
one of the eight possible (111) axes (expressed in the 
coordinates of the closely related cubic structure), with 
neighboring octahedra rotating in opposite directions in 
a cell-doubling, alternating three-dimensional 2x2x2 
checkerboard pattern. The crystal itself responds to the 
presence of these trigonal rotations about a common axis 
by stretching or compressing along that axis, forming a 
rhombohedral Bravais lattice. Selection of a specific oxy- 
gen octahedron as reference then permits eight possible 
distinct bulk configurations, each characterized and enu- 
merated by the particular choice of one of the eight pos- 
sible trigonal rotation axes for that specific octahedron. 
These eight reconstructions can be characterized either 
as ±[±1,±1, 1], namely a selection of an overall ± sign 
and a choice of one of the four unsigned [±1, ±1, 1 ] rota- 
tion axes, or, alternatively, simply as one of eight signed 
[±1,±1,±1] axes. 

The former perspective is visualized in Figure [3] The 
four (unsigned) [±1, ±1, 1 ] rotation axes designate the re- 
construction at an arbitrary origin of the crystal, where 
a final choice of sign determines whether the rotation 
at the origin is either clockwise ( "positive" ) or counter- 
clockwise ( "negative" ) , thereby fully specifying the mi- 
croscopic state of the crystal. Note that a change in this 
sign choice corresponds precisely to a rigid translation of 
the crystal by a [111] primitive translation vector (to a 
position in the crystal with opposite sign of rotation, due 
to the checkerboard pattern of the reconstruction) . Such 
a translation would not be observable on a macroscopic 
level. 
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FIG. 4: Sign pattern of octahedral rotations on opposite sides 
of stacking fault: octahedra in lower layer 21 (O) and upper 
layer 25 (□) and standardized displacement D ( — Pan- 
els (a) and (b) depict the two distinct choices for sign of ro- 
tation (+ and — , respectively) of 0<b ■ 



FIG. 2: Perovskite reconstruction in a~ a~ a~ Glazer notation, 
showing rotations of oxygen octahedra in opposite directions 
in a cell-doubling, alternating three-dimensional 2x2x2 
checkerboard pattern. (Degree of rotation exaggerated for 
illustrative purposes.) 




FIG. 3: "Positive" senses of the four (unsigned) [±1,±1,1] 
rotation axes for oxygen octahedra (§3 axis out of the page). 



A. Isolated stacking faults 

The RP phases observed in experiment d 49 * 52 *^) consist 
of superlattices of bulk perovskite slabs separated by the 
insertion of excess (001) AO planes to create stacking 
faults of the form . . ./AO/B0 2 /AO/AO/B0 2 /AO/. . . in 
the normal . . ./AO/B0 2 /AO/B0 2 /AO/. . . stacking se- 
quence. Across each such stacking fault, the bulk per- 
ovskite slabs are alternately displaced by in-plane vectors 
of the form ^f[±l, ±1, 0]. To our knowledge, however, no 
one has yet explored the effects of different combinations 
of possible symmetry-related bulk reconstructions on op- 
posite sides of the repeated AO planes. 

To enumerate the distinct possible configurations for 
such stacking faults, we restrict our consideration to cases 



where the material on cither side of the stacking fault 
possesses a specific bulk reconstruction throughout. We 
first focus on the reconstruction of the bulk material on 
the side "below" (at lower values for £3) the AO/AO 
stacking fault, which we shall denote as side 21. As dis- 
cussed above, the microscopic configuration of bulk ma- 
terial on side 21 can be fully specified by noting the ro- 
tation axis of the octahedron at an arbitrary, but from 
then onward fixed, origin 0%. The rotation axis of this 
particular reference oxygen octahedron can then assume 
any of eight choices among the (111) axes, which we can 
regard as a selection of an overall ± sign ("positive" or 
"negative") and a choice of one of the four (unsigned) 
[±1,±1,1] rotation axes. This specification then deter- 
mines the entire microscopic structure of the bulk mate- 
rial in side 21, according to the 2x2x2 three-dimensional 
checkerboard pattern of alternating signs for the rotation 
axes. For the present purpose, we choose always to se- 
lect the reference octahedron 0<& from among those octa- 
hedra immediately neighboring the stacking fault which 
have a positive rotation in the above sense of the choice 
of overall ± sign. 

Next, as noted above, the bulk perovskite on the other 
side, 03, of the stacking fault is displaced, in general, 
by a vector ^[±1,±1,0] relative to side 21. Temporar- 
ily disregarding the rotational state of all octahedra, we 
note that all four of these displacements are equivalent, 
since the octahedra on 03 are positioned at the centers of 
squares formed by the octahedra on 21. (See Figure [1]) 
This equivalence (apart from rotational states of the oc- 
tahedra) allows us to choose a standardized displacement 
D = ^[+1,+1,Q] from O a to select the origin O s . 
(More specifically, 0-% — 0% + T) + (a e 3 , where C ~ 3/2 
represents the vertical displacement between octahedra 
on opposite sides of the stacking fault, and S3 is a unit 
vector in the [001] direction.) 

The rotation state of 0<s then completely specifies the 
entire microscopic structure of the bulk material in side 
03. This rotation state can be specified as a selection of 
an overall ± sign ( "positive" or "negative" ) and a choice 
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of one of the four (unsigned) [±1,±1,1] rotation axes. 
Figure [4] illustrates exactly these two possible choices of 
sign, where the reference octahedron on side 03 is either 
positive (Figure a)) or negative (Figure [4jb>) ) . Generi- 
cally below, we shall refer to these two possible sign pat- 
terns in the stacking-fault configuration as "+" and "— " , 
respectively, reflecting the overall sign of the rotation of 
Oig. The combination of these two sign patterns and 
the four possible unsigned rotation axes (not indicated 
in Figure |1| for each origin, Oa and 0<b , leads to a total 
of 2 x (4 x 4) = 32 distinct possible configurations for 
this stacking fault. 

Next, we consider equivalence of stacking faults under 
application of C^ z symmetry. As discussed above, the ro- 
tation axis of O21 can always be selected to be among the 
four (unsigned) [±1,±1,1] axes. Because C± z symme- 
tries interconvert all four of these axes (C4 Z o[lll] = [111], 
Cf z o[lll] = [111], C| 2 o[lll] = [111]), an overall rotation 
of the coordinate system can be found to make the O53 
rotation about the [111] axis. We may thus define the ro- 
tation axis of 0% to always be the (unsigned) [+1, +1, 1 ] 
axis. This then narrows the phase space of thirty- two 
configurations listed above to now only 2x(lx4) = 8 
distinct configurations. 

These eight distinct configurations can be enumerated 
using two related nomenclatures (seen in Table [I]) , one 
which is algebraically explicit and suitable for symme- 
try arguments and the other which is more compact and 
convenient for communication. In the former case of the 
algebraic specifier, we enumerate each configuration by 
specification of the rotation states of each reference oc- 
tahedron, 0% and Org 1 expressing the rotation states in 
terms of (signed) [±1,±1,±1] vectors. For example, if 
O a is in rotation state [111] and is in rotation state 
[111], we write (111 ==> 111). For a more compact no- 
tation, we can take, without loss of generality as shown 
above, the rotation state of Oa to always be [111]. The 
unsigned rotation axis of 0<s is then related to that of 
0% by one of the four C*4 Z rotations of angles {0, J, ir, 
^}. To denote the configuration of the stacking fault, we 
then append the remaining choice of ± sign for 0<s to its 
rotation angle to produce a compound rotation-sign spec- 
ifier. For example, the (111 ==> 111) configuration from 
the preceding example may be more compactly written 
as [tt-] (since -C\ z o [111] = -[III] = [111]). Finally, 
Table [I] enumerates all eight configurations according to 
both their rotation-sign and algebraic specifiers. 

The above enumerated eight stacking-fault configura- 
tions can be reduced yet further by considering the ef- 
fects of additional, more subtle, symmetry operations — 
specifically, the Cix, C2y, and C$ z rotations and inver- 
sion /. To demonstrate the effects of such symmetry 
operations, we first consider two of the above eight stack- 
ing-fault configurations, [f + ] in Figure^ a) and [|-~] in 
Figure [5jb). From the figure, these two configurations 
clearly differ only in the overall choice of sign pattern for 
the upper layer 03 (n's). We now use symmetry argu- 
ments to demonstrate that these two stacking-fault con- 



Rotation-sign Specifier Algebraic Specifier 

[0+] (111 111) 

[0-] (111 => 111) 

[f+] (111 => 111) 

[§"] (111^111) 

[tt+] (111 => 111) 

[TV ] (111 111) 

[4f+] (111 =► 111) 

[f-] (111 => 111) 



TABLE I: Enumeration of interfacial stacking faults remain- 
ing after application of C*4 Z and translational symmetry, la- 
beled by either relative rotation and sign of the bulk regions 
on opposite sides of the stacking fault (rotation-sign specifier) 
or rotation state of the two reference octahedra (algebraic 
specifier). 




FIG. 5: Full specification of octahedral rotations on opposite 
sides of stacking fault: octahedra in lower layer 21 (O) and up- 
per layer 25 (□) and standardized displacement D ( — ► ), with 
explicit (signed) [±1, ±1, ±1] rotation axes to specify full con- 
figuration. Panels (a) and (b) depict different configurations, 
with Panel (c) showing the application of Ci x rotation to (a) 
to produce a configuration equivalent to (b). 

figurations are actually equivalent, related simply by a 
C2X rotation. 

The general strategy which we shall employ is to ex- 
ploit the fact that the algebraic specifier not only deter- 
mines the rotation state of all of the octahedra in the 
crystal but also transforms in relatively simple ways un- 
der the application of symmetry operations to the overall 
crystal. Figure [5^c) illustrates the state of the crystal 
in Figure [5^ a) after the application of a C<i x symmetry 
about a specific axis. To ensure that the basic crystalline 
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structure maps back onto itself, this axis is chosen to 
pass at equal distances to the two neighboring AO layers 
in the stacking fault and through the point immediately 
above the 0% reference octahedron (base of the standard- 
ized displacement D in Figure [5ja,) ) . The application of 
Ci x effects the following changes to the configuration in 
Figure [BJa) : 

1. the position of 0% maps from the bottom layer 21 
to the top layer 25 (this changes the central O in 
Figure [5ja) into the central □ in Figure |5jc)); 

2. the rotation state of O3 maps from [111] to [111]; 

3. the position of the 25-side octahedron (□) in the 
lower-right quadrant of Figure |5ja) maps to the 21- 
side octahedron (o) in the upper-right quadrant of 
Figure ^c); 

4. the rotation state of this latter octahedron maps 
from [111] to [111]. 

To determine the algebraic specifier for the final config- 
uration in Figure |5jc), we now identify an octahedron 
on the lower side of the stacking fault which has a posi- 
tive rotation from among the (unsigned) [±1, ±1, 1 ] axes 
(base of the standardized displacement D in Figure[5jc)). 
The rotation state of this octahedron, combined with 
that of the corresponding reference octahedron on the 
upper side, defines the algebraic specifier for the result- 
ing stacking fault. For this particular case, it is evident 
from the figure that C 2x o [f +] = (111 => 111). Ref- 
erence to Table [j] then identifies (111 =S> 111) = [§"], 
so that C 2x [f ] = [§ _ ], as direct comparison of Fig- 
ures ^Jo) and (c) confirms. Thus, we conclude that [^ + ] 
and [§~] are identical and related by the C 2x rotation. 

To develop an algebra for the action of symmetry oper- 
ations on arbitrary stacking-fault configurations, we con- 
sider the effects of such operations on general stacking- 
fault configurations, (abc => afl-y). For the case of Ci x , 
we note in preparation that the action of Ci x on any 
vector (xyz) is C2X (xyz) = (xyz). Then, as above, the 
operation of C2X transforms the generic algebraic speci- 
fier (abc a/?7) as follows: 

1. it swaps abc and oePj; 

2. it maps abc to abc; 

3. it changes the sign of a/37; 

4. it maps a/?7 to a/3"f 
In short, 

Cix (abc => ctft^f) = (a/?7 abc). 

Table [TT| summarizes the effects, similarly determined, on 
stacking-fault configurations of a number of useful sym- 
metry operations: rotations C 2x , Qy, C 2z , and C^ z and 
inversion /. (Note that inversion, through the midpoint 



Algebraic Specifier 
Symmetry 



Original Final 

C2x (abc => a/37) ( a Pl a bc ) 

C2 V (abc =>■ a/37) ( a Pl =^ abc ) 

C2z (abc => a/37) ( abc => a/37) 

Caz (abc => a/37) ( bac f3a^) 

I (abc => a/37) ( a Pl =^ abc ) 



TABLE II: Rotation (C) and inversion (I) symmetry opera- 
tions applied to general stacking-fault configuration. 

of the vector connecting the two reference octahedra Oqi 
and 0<s, requires special care: it does not change the 
direction of pseudo-vectors such as the rotation axes of 
the different octahedra but simply interchanges the two 
sides, 21 and 25.) 

We are now able to demonstrate further symmetry re- 
duction of the stacking-fault configurations. Table [I] lists 
all eight possible stacking fault configurations, provid- 
ing both rotation-sign and algebraic specifiers for each. 
Application of the symmetry operations from Table [TT] 
to these configurations, and subsequent translation T by 
ao[100] if necessary to make Oa of rotation state [111], 
generates the following relationships, 

ToC 2a o[f+] = [f-] (la) 
T°C 2x o[f +] = [§-] (lb) 
C 4z o/o[f+] = [f-]. (lc) 

The above symmetry operations, which each involve in- 
terchange of the 21 and 25 sides of the interface, demon- 
strate equivalence among the set of four stacking faults, 
{[f + L [f~L [ir + ]> [lT~]}- Thus > only five unique stack- 
ing-fault configurations remain, [0 + ], [0 - ], [f + ], [tt + ], 
and Comparison of the algebraic specifiers of these 

five configurations, from Table [T] above, with the com- 
puter-generated distorted phases listed as Nos. 18-22, in 
Table III of Hatch et al.pl confirms that our result is con- 
sistent with theirs, though theirs is limited to the n = 1 
RP phase. We confine ourselves to discussion of these 
five unique configurations for remainder of this work. 



B. RP phases 

The RP phases, which we study in this work, consist 
of periodic arrays of the above configurations of stack- 
ing faults. While constructing such arrays from individ- 
ual stacking faults, the array periodicity may be main- 
tained either through simple alternating patterns of bulk 
regions (21/23/21/25 sequencing) or through more com- 
plex, and possibly lower-energy, patterns (for example, 
2l/2ye:/2l/2V£ sequencing). 

Immediately below, we demonstrate that, in fact, an 
21/23/21/25 sequence of bulk regions always corresponds 



8 



to repetition of the same type of stacking-fault config- 
uration. Three very compelling reasons then follow for 
studying RP phases with this type of periodicity. First, 
one should expect the preferred RP phase to consist of a 
sequence of stacking faults, which are all of the lowest-en- 
ergy configuration; this would then naturally correspond 
to an 21/05/21/58 sequence. Second, actual experiments 
find super lattice sizes consistent with small, simple re- 
peat units.^- 52 * " lFinally, isolated stacking faults are also 
observed under certain conditions,^ and the study of RP 
phases of 21/05/21/05 periodicity, containing two identi- 
cal faults in each primitive unit cell, then allows for the 
extraction of the behavior of individual isolated faults. 

To establish that an 21/58/21/58 sequence of bulk re- 
gions corresponds to repetition of the same configuration 
of stacking fault, we must prove the equivalence of the 
two stacking faults, ( 21/05) and (05/21)', in such an RP 
phase. As in Section |III A| we take the algebraic speci- 
fier for this first stacking fault (21/58) to be the generic 
(abc =>• ctpj) where abc and a/?7 refer to the rotation 
states of the two reference octahedra for this fault, Oa 
and Org, respectively. We must then choose two corre- 
sponding reference octahedra and for the second 
stacking fault (58/21)'. Since D = D, we can set = 
0% + (n — l)£aoe3 and = 0<g, + (n — l)(ao<§3, where n 
is the number of the RP phase (A n+ iB n 03 n +i) so that 
n— 1 layers of material separate Oa from and similarly 
separate Org from . (Consideration of the n — 1 case 
should make this apparent.) The alternating octahedral 
rotations of this a~a~a~ Glazer system thereby define 
the algebraic specifier for the (58/21)' stacking fault as 
((-l) n - x (aj97) =*> (-ly^iabc)). However, displace- 
ment of the arbitrary origin (reference octahedra) by 
ao[100], for cases when n is even, generates an equivalent 
algebraic specifier for this (58/21)' interface which is now 
insensitive to n, (a[3"f abc). Finally, we apply inver- 
sion symmetry, as introduced in Table pT| to this (58/21)' 
stacking fault, I o (a/37 == ^ > a ^ c ) = (abc =>• a(3j), and 
thus prove that the (21/58) stacking fault is equivalent to 
the (05/21)' fault. 

Lastly, since all stacking faults in an 21/05/21/05 RP 
phase always possess equivalent configurations, we can 
uniquely label these RP phases using the same rotation- 
sign specifiers established above for the five symmetry- 
distinct stacking faults: [0 + ], [0~], [f ], br + ], and [ir~]. 



IV. MODEL AND COMPUTATIONAL 
METHODS 

A. Shell potential 

As discussed in Section |TT] above, to study the generic 
behavior of a~a~a~ perovskites which form RP phases, 
we employ a shell- pot ential modeP^ parameterized for 
strontium titanateP^ Shell-potential models are formu- 
lated as an extension to ionic pair potentials and em- 
ployed to capture the polarizability of the atomic con- 



stituents. The shell model separates each ion into two 
parts, a core and an outer shell, which possess individual 
charges that sum to the nominal charge of the ion. The 
total model potential U consists of three terms, 



U= Up + Uc + Ub, 



(2) 



representing, respectively, the polarizability of the ions, 
and the Coulomb and short-range interactions among the 
ions. The polarizability is captured by harmonic springs 
connecting the core and shell of each ion, so that Up has 
the form. 



U P = 



(3) 



where |Arj| is the core-shell separation for ion i and the 
ki are a set of ion-specific spring constants. The Coulomb 
contributions take the form, 



Ur 



' k c q,qj 



(4) 



where i and j range over all cores and shells (excluding 
terms where i and j refer to the same ion) , and qj are 
the corresponding charg the distance between the 

charge centers, and k c is Coulomb's constant. Finally, 
the short-rang e int eractions are included through a sum 
of Buckingham 1 -^ pair potentials (which can be vie wed 
as combinations of Born-Mayer^nSl and Lennard- Jones^H 
potentials) of the form, 



A, 



i/Pi. 



C ij r ij 



(5) 



where i and j range over all shells and A+j, pij, and CV, 
are pair-specific adjustable parameters. Here, the first 
term (Born-Mayer) serves as a repulsive short-range in- 
teraction to respect the Pauli exclusion principle, and 
the second term (Lennard- Jon es) m odels the dispersion 
or van der Waals interactions P^Sl The specific electro- 
static and short-range shell-model parameters used in 
thi s st udy were fit to strontium titanate b y A khtar et 
a j Jiou with values as listed in Tables III and IV Finally, 



we wish to emphasize again, as it is rarely mentioned 
explicitly in the shell-potential literature, that the pair- 
potential terms in Ub apply to the shells only, and not 
to the cores. 

Shell models have been extensively used for decades as 
the primary empirical potential for modeling perovskites 
and other oxidesP ^ 6 * 107 ! We tested the correctness of our 
coded implementation of this potential through compar- 
isons of lattice constants and elastic moduli without the 
AFD reconstruction and find excellent agreement. For 
instance, using the same shell potential and the same 
non-reconstructed ground state, we predict a volume per 
SrTiC-3 chemical unit of 59.18 A 3 , which is within 0.4% 
of the value calculated by Akhtar et alP^H For the elastic 
moduli, we find Cu = 306.9 GPa, C 12 = 138.7 GPa, 
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Shell 


Core 


Spring Constant 


Ion 


Charge [e] 


Charge [e] 


[eV-A- 2 ] 


Sr 2 + 


1.526 


0.474 


11.406 


Ti 4+ 


-35.863 


39.863 


65974.0 


o 2 - 


-2.389 


0.389 


18.41 



TABLE III: Electrostatic shell-mod el p otential parameters 
used in this study (from Akhtar et alpSU). 



Interaction A [eV] p [A] C [eV ■ A 6 ] 

Sr 2+ ^ O 2 - 776.84 0.35867 0.0 

Ti 4+ <^ O 2 - 877.20 0.38096 9.0 

O 2 " <^ 2 ~ 22764.3 0.1490 43.0 



TABLE IV: Short-range shell-mode l po tential parameters 
used in this study (from Akhtar et alP^H). 

and C 44 = 138.8 GPa, which are within 1.8%, 1.0%, 
a nd .7%, respectively, of the values from Akhtar et 
alp 3 From this, we conclude that our implementation 
of the potential is correct. We also would like to note 
that, when the a~a~a~ reconstruction is considered, 
significant changes occur, and we find, instead, a vol- 
ume per SrTiC>3 chemical unit of 58.49 A 3 and elastic 
moduli of Cxi = 275.6 GPa, C 12 = 144.3 GPa, and 
C44 = 133.5 GPa. This underscores the importance of 
considering AFD reconstructions when constructing such 
potentials. 

B. Numerical methods 

In this work, we compute the Coulombic interac tion^^ 
from Q using a Particle Mesh Ewald algorithm^SHlliJ 
with all real-space pair-potential terms computed out to 
a fixed cutoff distance using neighbor tables. Analytic 
derivatives are used to determine the forces on the cores 
and shells within the supercell, and finite differences are 
used to compute the generalized forces on the superlat- 
tice vectors. While such finite-difference methods occa- 
sionally introduce issues with numerical precision into 
the resulting calculations, we are able to mitigate such 
effects with care in selection of the finite-difference step, 
scaling it proportionally to the sizes of the lattice vectors 
involved. 

With the gradients determined as above, we relax 
each system using the techn ique of preconditioned con- 
jugate gradient minimization (specifically, the Polak- 
Ribier d 113 ! method), fully optimizing the ionic coordi- 
nates in an inner loop of the routine and then optimizing 
the lattice vectors in an outer loop. The preconditioner 
applies only to the ionic relaxation of this procedure and 
scales the generalized force on the geometric center of 
each ion (mean position of the core and shell) separately 
from the generalized force on the core-shell displacement, 
with the latter scaled in inverse proportion to the spring 



constant of the corresponding ion. This process ensures 
that all supercell energies are relaxed with respect to both 
ionic and lattice coordinates to a precision of ~0.3 u.eV 
per supercell. 



C. Ground state 

To establish the ground state of the model potential 
used in our calculations, we have carried out what we 
regard as a thorough, but not exhaustive, search for a 
probable ground-state structure. Indeed, we have found 
no alternative structure which relaxes to an energy less 
than our candidate ground-state structure within our po- 
tential. We performed quenches on hundreds of random 
displacements from the idealized positions of the lxlxl 
primitive unit cell to explore various potential reconstruc- 
tions for supercells up to 6 x 6 x 6. We also considered a 
number of highly ordered configurations commensurate 
with antiferrodistortive disordering. 

Among those minima which we explored, we selected 
the lowest-energy configuration to serve as the bulk crys- 
talline state throughout this study. This configuration 
possesses a fairly regular pattern as depicted in Fig- 
ure [2j namely each oxygen octahedron rotates slightly 
along trigonal directions in a cell-doubling, alternating 
three-dimensional 2x2x2 checkerboard pattern, which 
is precisely the desired a~a~a~ Glazer system from Sec- 
tion HU We also find that the lattice vectors assume the 
rhombohedral symmetry commonly found in a~a~a~ 
perovskites, with the cubic lattice stretching along the 
same trigonal axes about which the octahedra rotate. 

The specific lattice vectors which we find for the 2 x 2 x 
2 reconstructed supercell of our model are the columns 
of the following matrix, 

/ 7.764 0.028 0.028 \ 
R = 0.028 7.764 0.028 A. 
\ 0.028 0.028 7.764 / 

D. Construction of RP supercells 

Our construction of initial ionic configurations for the 
A„+iB n 3 „ + i RP phases proceeds through a detailed 
process defined in Section |III| As discussed there, we 
carefully consider a large number of combinatorial pos- 
sibilities resulting from the different geometries due to 
the aforementioned 2x2x2 cell-doubling, a~a~a~ bulk 
reconstruction. We construct each configuration of the 
RP phase by stacking 2n layers of bulk perovskite along 
the [001] direction, with the two stacking faults (excess 
AO (001) planes) so situated as to equally divide the 
supercell into two bulk slabs, each with n layers and po- 
tentially different a~a~a~ reconstructions. Relative to 
each other, these bulk slabs are also laterally displaced 
by one of the four (equivalent) ^[±1,±1,0] vectors at 
each stacking fault. We perform this procedure to create 
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all eight (111 af3j) configurations from Section III A 



through combinations of reconstructions in the two bulk 
slabs, for RP phases from n = 1 ... 30, containing from 56 
to 1216 ions, respectively. Each RP superlattice system 
is then fully relaxed, as described above, both in terms 
of internal ionic coordinates and lattice vectors, until the 
energy is minimized to within a precision of ~0.3 u.eV. 
In all cases, we find that the material on either side of 
the stacking fault maintains a uniform bulk reconstruc- 
tion throughout, so that the enumeration of Section [ill A| 
is suitably complete. Finally, we test for meta-stability 
by creating ensembles of 20-250 samples for each config- 
uration of octahedra rotations for all RP phases up to 
n = 10, introducing root-mean-square displacements of 
0.001 A for each coordinate axis of every core and shell, 
and quenching each of these randomized systems. 



V. RESULTS 

As described in Section |TT] above, the results obtained 
below have been computed with the standard shell-po- 
tential model for strontium titanate ) 101 * 106 * 10 ^ but are 
interpreted to represent the generic behavior of an anti- 
ferrodistortive perovskite (ABO3) from among the class 
of perovskites within the a~a~a~ Glazer system which 
form RP phases. 



A. Formation of RP phases 

After constructing the various possible RP phases ac- 
cording to the prescription in Section [IVD| the ideal pro- 
cedure for determining the energetics of the RP phases 
is to fully relax each system, in terms of both internal 
atomic coordinates and superlattice vectors, without any 
externally imposed symmetries. The energies of the ref- 
erence bulk materials, perovskite (ABO3) and A-oxide 
(AO), should also be fully relaxed in the same sense. 

Without full relaxation of lattice vectors for both the 
RP phases and bulk reference cells, any attempts to ex- 
tract energies will introduce significant errors that scale 
in proportion to the amount of bulk material between 
stacking faults. (Such failures to fully relax the lattice 
vectors might arise due to expediency in ab initio stud- 
ies or through attempts to simulate higher-temperature 
phases by imposing a higher-temperature lattice struc- 
ture while relaxing ionic coordinates.) 

Figure |6^a) shows actual results obtained from such 
a methodologically flawed procedure, where the energy 
needed to form the RP phase from the bulk perovskite 
ABO3 and bulk A-oxide AO, Sa„ + iB„o 3 „ + i ~ n ■ Eabo 3 — 
Ej^q, is computed without full relaxation of the lat- 
tice vectors of the RP phase, but instead with the en- 
forcement of tetragonal lattice vectors associated with 
a higher-temperature symmetry not possessed by the 
ground state of our model. Under these constraints, the 
energy exhibits a clear linear increase with n, the amount 
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FIG. 6: Formation energies of RP phases per excess AO chem- 
ical unit as a function of n for various members of homolo- 
gous series A n +iB n 03 n +i: methodological error of employ- 
ing (a) tetragonal supercells for RP phases or (b) fully re- 
laxed supercells for RP phases but unrelaxed lattice vectors 
for bulk reference cells. Each panel shows results for the dis- 
tinct stacking-fault configurations corresponding to different 
arrangements of octahedral rotations, as enumerated in the 
text (Section [UTA I. 



of bulk material between stacking faults, due to an in- 
flated RP-phase energy associated with artificial strain 
in the bulk regions. (Note that the results in the figure 
are visually indistinguishable whether or not one allows 
the perovskite or A-oxide bulk reference cells to fully re- 
lax.) Figure[6](b) shows results obtained where the lattice 
vectors of the RP phases are allowed to relax fully but 
the lattice vectors of the bulk reference cells are not. In 
this case, excess bulk energy is not included in the RP 
cell but instead in the bulk reference cells, and the en- 
ergy of formation now decreases linearly with n. Either 
type of error will confound the extraction of meaningful 
information on the RP phases. The calculation of such 
phases is herein seen to be quite delicate, a fact which 
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FIG. 7: Formation energies of RP phases per excess AO chem- 
ical unit as a function of n for various members of homolo- 
gous series A n +iB„03 n +i: labeled curves represent each of 
the stable stacking-fault configurations corresponding to dif- 
ferent arrangements of octahedral rotations, as enumerated 
in the text (Section III A I. Results calculated using fully re- 
laxed ionic and lattice coordinates for both RP phases and 
bulk references. 



only becomes obvious when calculations are pursued for 
large values of n. 

Figure [7] presents the results, properly calculated with 
full relaxation of all lattice vectors, for the energies of all 
eight (111 a/37) stacking-fault configurations of our 
model potential, relative to the corresponding amounts of 
the fully relaxed, reconstructed perovskite and the fully 
relaxed A-oxidc material. As in the examples above, 
the energies are expressed per excess AO chemical unit, 
2 x 2 = 4 of which appear in each reconstructed RP su- 
perlattice unit cell (and 2 x (2 x 2) = 8 of which appear 
in each supercell with two stacking faults and thus two 
RP super lattice unit cells). The final relaxed RP con- 
figurations indeed are not distinct, as expected from the 
symmetry arguments in Section [ill A[ which showed that 
[f +] = [f "] = [^ + ] = [^f "], so that there should be at 
most five different configurations. Somewhat unexpect- 
edly, rather than displaying five separate curves for the 
five symmetry-distinct configurations ([0 + ], [CP], [§ + ], 
[tt + ], [7r~]), Figure [7] instead shows only four curves: we 
find that the [| + ] configuration is unstable and relaxes 
directly to the [tt~] configuration, as determined by en- 
ergy calculations and observations of transitions in lattice 
symmetries for each configuration during full relaxation. 
Additionally, the [CP] configuration appears meta-stable, 
often tran sitioni ng, upon small perturbations described 
in Section IV D to the [0 + ] configuration. The remain- 



der of this discussion thus only reports on the behavior of 
the final three distinct, stable configurations ([0 + ], [tt + ], 
[tt~]) as well as the one additional distinct, meta-stable 
configuration ([0 - ]). 

From Figure [7] we note immediately that, depending 
upon the relative rotational orientation of oxygen octa- 



hedra on cither side of the stacking faults, the config- 
urations separate into two groupings, one with positive 
energy and one with negative energy relative to phase 
separation into bulk perovskite and bulk A-oxide. The 
energy scales associated with the different configurations 
may thus be sufficient to either stabilize or destabilize 
the formation of an RP phase. 

Moreover, for our model potential, the two more stable 
stacking-fault configurations ([0 + ] and [tt~]) possess neg- 
ative energies for all n, providing the first demonstration 
that RP phases can be favored over formation of bulk 
perovskite and A-oxide for all values of n. Also, we find 
the [tt~] configuration to be favored over all other pos- 
sible configurations in our model, regardless of the value 
of n. 

It is also immediately apparent that the RP phases for 
n = 1 deserve special consideration for our model poten- 
tial. The n — 1 phases for the two configurations with 
positive energy collapse to related n = 1 configurations 
with negative energy, with [tt + ] becoming [tt~] and [CP] 
becoming [0 + ]. This might be related to the realization 
that the full 2x2x2 reconstruction cannot be expressed 
on either side of the stacking fault for an n = 1 RP phase, 
and thus these higher-energy configurations forfeit some 
stability for the n = 1 RP phase. 

The magnitude of the energy differences between con- 
figurations is not only sufficient to affect the stability of 
the phases, but is also significant compared to thermal 
energy scales. For instance, we find the separation be- 
tween the stable and unstable groupings to be ~200 meV 
per excess AO chemical unit (or ^800 meV per RP su- 
perlattice unit cell, (2x2)- A n+1 B„03„+i). To place this 
energy in context, we can also consider the energy per in- 
terfacial oxygen octahedron (one on the top and bottom 
of the stacking fault for each excess AO chemical unit), 
^100 meV, which is about ~30 meV per configurational 
(rotational) degree of freedom. This is larger than, but 
comparable to, the thermal energy at room temperature, 
26 meV, so that at significant temperatures we still can 
expect a small, but statistically noticeable, tendency for 
the stacking fault to preferentially select the minimal en- 
ergy configuration over other disfavored configurations. 

This energy scale is also significantly larger than the 
energy scales we might expect for excitations in the bulk 
of the material. For instance, the structural phase transi- 
tion in strontium titanate (the mate rial to which the orig- 
inal model actually had been fitP^U) occurs near 105 Kpl 
which corresponds to an energy scale of 9 meV. To fur- 
ther place the above ~30 meV interfacial energy scale 
in context, we consider the energies within our model 
of two alternate rotational 2x2x2 bulk reconstruc- 
tions, those with rotations about (110) axes (orthorhom- 
bic) and (100) axes (tetragonal), finding rotational en- 
ergy scales of 4 meV and 16 meV, respectively, relative to 
the ground-state reconstruction. These comparisons lead 
to the intriguing conjecture that, for some perovskites 
with low densities of stacking faults (high n) , there exist 
a range of temperatures above the bulk structural phase 
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transitions, in which the bulk rotational energy scale is 
insufficient to constrain the rotations of the octahedra, 
while the interfacial energy scale is suitable to constrain 
the octahedral rotations to specific preferred orientations. 
This corresponds to a picture of "fuzzy" octahedra in the 
bulk regions with noticeable orientational preferences re- 
maining for the octahedra at the stacking faults. Alter- 
natively, for perovskites in RP phases with low values 
of n and higher densities of stacking faults, one might 
expect this increased energy scale to raise the transition 
temperature associated with octahedral rotations. 



B. Stacking fault interactions and 
intergrowth formation 

While Section |V A| discussed the relationship among 
energies of formation of different configurations, it is fur- 
ther instructive to examine the energy and behavior of 
each configuration independently, as a function of sep- 
aration between stacking faults. Figure |8] plots the in- 
teraction energies per excess AO chemical unit for each 
configuration, which, for a given configuration, are the 
formation energies of each phase relative to the energy 
of its n = 2 phase. (For the two stable, low-energy con- 
figurations within our model, the n = 2 phase is, in fact, 
lowest in energy, so that this plot represents the binding 
energies of the stacking faults for these configurations. 
For the two high-energy configurations, unstable to dis- 
sociation into separate regimes of bulk perovskite and 
A-oxide, the n = 1 phase is actually lowest in energy, 
so that this plot approximately represents the binding 
energies of the stacking faults.) For our model, in all 
cases, there is a relatively weak, but consistently attrac- 
tive, interaction (apart from special cases for n = 1), 
varying by over a factor of two from ^7-14 meV per ex- 
cess AO chemical unit, depending upon the stacking-fault 
configuration. The strengths of the interaction occur in 
precisely the same sequence as the stability of the config- 
urations, with the more attractive interactions occurring 
in the more stable configurations. 

Within the current model, the universally attractive 
and eventually asymptotic nature of the interaction for 
n > 2 leads to a set of interaction energies with down- 
ward curvature. The standard secant construction for 
mixed phases then indicates that, for n > 2, the en- 
ergy of a system can always be lowered by forming a so- 
called "intergrowth" phase, mixing regions of lower- and 
higher-n phases (i.e., non-uniformly changing the sepa- 
ration between stacking faults that define the phases) 
while maintaining a consistent average stoichiometry. To 
understand how this arises from the particular form of 
these interactions, note that changing the location of a 
particular fault (in a uniform n > 2 phase) decreases the 
distance to one of its neighboring faults but increases the 
distance to the other neighboring fault. The downward 
curvature of the interaction energy, however, ensures that 
the energy gain from the decreased distance more than 
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FIG. 8: Interaction energies among stacking faults (forma- 
tion energies relative to n — 2 RP phase) per excess AO 
chemical unit as a function of n for the homologous series 
A n +iB n 03n+i (same notations as in Figure [7b. Data depict 
relatively weak, but consistently attractive, interactions in all 
cases. 



compensates for the cost from the increased distance, so 
that the overall energy is lowered by varying the distances 
between faults. 

Finally, we consider the possibility of intergrowth phe- 
nomena. The two higher-energy configurations within 
our model are unstable with respect to dissociation into 
inter growths of bulk perovskite and bulk A-oxide, and 
thus need not be considered further here. To examine the 
possibility of intergrowths in the two lower-energy config- 
urations, we consider environments both with a relatively 
low density (high n) of stacking faults and those with a 
higher density (lower n) of stacking faults. In the former 
case, the stacking faults attract to form dense phases at 
the value of n for which the energy is minimized, namely 
ri = 2 within our model, alternating with regions of bulk 
perovskite. In the latter case, the additional A-oxide ex- 
ceeds the concentration required for creation of a consis- 
tent phase at the preferred minimum energy n = 2 value. 
For the stoichiometry such that 1 < n < 2 on average 
throughout our model system, the excess A-oxide could 
either form intergrowths consisting of the n = 2 phase 
and bulk A-oxide or intergrowths consisting of n = 2 and 
n = 1 phases. The creation of n = 1 phases, by insertion 
of excess AO planes, corresponds to the transformation 
of a given unit of the n = 2 phase into two units of the 
n = 1 phase. Since twice the formation energy of the 
n = 1 phase is always lower than the formation energy 
of the n = 2 phase within our model, we recognize that, 
for stoichiometries with 1 < n < 2 on average, the two 
stable configurations ([0 + ] and [ir~]) will always form in- 
tergrowths of n = 1 and n = 2 phases, rather than of 
bulk A-oxide and n = 2 phase. For RP phases with yet 
higher concentrations of species A (n < 1), a similar ar- 
gument shows that, within our model, intergrowths form 
between the n = 1 RP phase and bulk A-oxide. 
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FIG. 9: Interaction energies among stacking faults (forma- 
tion energies relative to n = 2 RP phase) per excess AO 
chemical unit as a function of 1/n for the homologous series 
A n +iB n 03n+i: extraction of energies to infinite superlattice 
size (dotted lines) to determine binding energies of stacking 
faults. Linear behavior indicates an inverse-distance interac- 
tion between stacking faults. 



Configuration 


Binding Energy [meV] 


[o-]t 


7.218 


[7T+] 


9.476 


[0+] 


13.288 


[TP"] 


14.658 


t Meta-stablo 



TABLE V: Binding energies of stacking faults per excess AO 
chemical unit for different configurations. 



C. Isolated stacking faults 

Although the stacking faults are attractive in our par- 
ticular model, crucial early growth experiments on RP 
phases often observed these stacking faults in isolation,^ 
perhaps due to kinematic constraints. Here, we consider 
the physical properties of such isolated faults. 

To extract the behavior of isolated stacking faults at 
infinite separation, we first determine the functional de- 
pendence of the interaction on the separation between 
stacking faults. Figure [9] shows the interaction energies 
of the stacking faults from Figure [8] now plotting these 
energies as a function of 1/n, where n measures the sep- 
aration between stacking faults in terms of the number 
of intervening bulk layers. The data in this figure exhibit 
clear linear behavior, indicating that the interaction is in- 
versely proportional to fault separation. The ordinate in- 
tercepts (1/n. — ► 0, or n — ► oo) denote the energy released 
per excess AO chemical unit when stacking faults bind 
together into RP phases at their most energetically fa- 
vored value of n (for the two stable configurations within 
our model, [0 + ] and with energy minima at n = 2). 
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FIG. 10: Formation energies of RP phases per excess AO 
chemical unit as a function of 1/n for the homologous series 
A n +iB n 03n+i: extraction of energies to infinite superlattice 
size (dotted lines) to determine formation energies of isolated 
stacking faults. 



Configuration 


Formation Energy [meV] 


[o-] J 


162.0 


[7T+] 


143.3 


[0+] 


-21.0 


[7T-] 


-44.5 



* Meta-stable 

TABLE VI: Formation energies of isolated stacking faults per 
excess AO chemical unit for different configurations. 



Table [V] summarizes the values of these binding en- 
ergies per excess AO chemical unit, as determined by 
linear least-square fits to the energies for RP phases in 
Figure [9j ignoring in each case the three smallest separa- 
tions (n < 3) which clearly involve significant higher-or- 
der interactions. Even for the most attractive configura- 
tion ([7r~]), our model shows a ~7 meV binding energy 
per octahedron neighboring the stacking fault (recall that 
there are two interfacial octahedra for each excess AO 
chemical unit, one on either side of the stacking fault), 
which is relatively low compared to room temperature, a 
fact which may be related to the observation of isolated 
faults in some materials. 

With the form of the interaction now determined, we 
turn to the formation energies of isolated stacking faults 
in their various configurations. We can now finally deter- 
mine whether these different configurations have a large 
effect on formation energies of RP phases, as Section |Y A| 
suggests. Figure [T0| shows the formation energy per ex- 
cess AO chemical unit as a function of 1/n. 

Here, the ordinate intercepts give the energy to in- 
corporate an AO chemical unit into an isolated, excess 
AO plane in perovskite relative to that chemical unit ap- 
pearing in bulk A-oxide. Table [VT| summarizes the for- 
mation energies of isolated stacking faults as determined 
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by the same type of linear least-square fit used to extract 
the binding energies above. Within our model, we again 
find significant differences in formation energies between 
the two groupings of configurations (those with positive 
and negative energies), of up to ^200 meV per excess 
AO chemical unit or ~100 meV per oxygen octahedron 
neighboring the stacking fault. We also identify a much 
smaller energy difference between the two lowest-energy 
configurations (those with negative energies) of 23.5 meV 
per excess AO chemical unit, or 11.8 meV per octahe- 
dron. Following a similar argument as in Section |V A| 
for a range of temperatures above the bulk phase tran- 
sitions, we may still find that octahedra near isolated 
stacking faults are constrained to a subset of low-energy 
rotational orientations, while those in the bulk fluctuate 
among all possible rotations. 

Finally, comparing the results in Table |VI| with those 
in Table [Vj we come to the commonsensical, though not 
logically necessary, result that the sequence of stacking 
faults from the most stable (i.e., lowest formation energy) 
to the least stable, directly corresponds to the sequence 
of binding strengths, from the strongest attraction to the 
weakest. We also propose a simple physical mechanism 
for differentiating between higher and lower energies of 
formation for the configurations, potentially reflecting 
their stability relative to phase separation into bulk per- 
ovskite and bulk A-oxide. Figure [TT] depicts the planar 
movement of oxygens in the AO layers on opposite sides 
of the stacking fault. The sense of rotation of the interfa- 
cial octahedra drives constituent oxygen ions to move in 
specific patterns within their respective AO layers along 
the stacking fault. For the high-energy configurations, 
[CP] or [tt + ], the interfacial oxygen ions from opposite 
AO layers are displaced toward each other, despite their 
natural Coulombic repulsion, as in F igu re [Til a). For the 
low-energy configurations in Figure 11 b), |0 + ] or 
the same oxygen ions move past each other along the 
interface, tending to lessen the approach of like-charged 
ions toward one another and thus minimizing the energy 
Thus, we find that the stability of different stacking-fault 
configurations can be readily understood in terms of sim- 
ple Coulombic effects. 



D. Point defects 

With the energies of the various RP phases now de- 
termined, we are finally in a position to substantiate the 
claim that, in A-rich stoichiometries of ABO3, the RP 
phases are preferred over the formation of point defects. 

Interstitials in perovskites are known to be energeti- 
cally very unfavorable.^ On the other hand, A-rich stoi- 
chiometries can be understood by noting that the ad- 
dition of an AO chemical unit is identical to the incor- 
poration of a bulk ABO3 chemical unit followed by the 
removal of a BO2 chemical unit, the latter step forming a 
BO2 vacancy complex. To decide whether this formation 
of vacancies is more or less favorable than the incorpo- 
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(b) Low-energy configurations: [0+] or [ir ] 

FIG. 11: Planar movement of oxygen ions on opposite sides of 
AO/ AO stacking faults: rotation of octahedra in lower layer 
21 (/y , outlined arrowheads on solid grid) and upper layer 25 
{/y 1 outlined arrowheads on dashed grid) projected onto in- 
terfacial plane and resulting movement of oxygens (*s\ , filled 
arrowheads) within such plane. Panels depict (a) high-energy 
configuration with movement of neighboring oxygen atoms 
directly toward each other and (b) low-energy configuration 
with oxygen atoms moving past each other. 



ration of excess AO chemical units into planar stacking 
faults of the type in RP phases, great care must be taken 
to accurately balance stoichiometries. 

In the current context, we have been considering the 
excess AO chemical units originating from an external 
reservoir of A-oxide. To maintain a consistent chemi- 
cal potential, the formation energy of an isolated BO2 
vacancy in bulk perovskite must be computed as the en- 
ergy first to remove an extra AO chemical unit from bulk 
A-oxide, and then to insert it into a large supercell con- 
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FIG. 12: Formation energies of BO2 vacancies as a function 
of inverse-square supercell dimension 1/rf: extraction of en- 
ergies to infinite supercell (dotted line) to determine energy 
of isolated vacancy. (See text for definition of 77.) 



taining n chemical units of ABO3. Expressed mathemat- 
ically, the energy of formation of this vacancy is 



Ev — —Eao + £-A„ + iB„0 3 „ + i — ^A„B„0 3 „ 

= —Eao + ^A m B m _i0 3m _2 — Ea„ 1 _ 1 b„ 

= £ ; A m B m _ 1 3m -2 - [( m - 1)-EaB0 3 • 



(6) 



-i«j 3 ( m 

E AO ] 



where m = n + 1. As the final equality in (|6| indicates, 
the energy of formation of an isolated BO2 vacancy in 
bulk perovskite is equivalent to the energy of a BO2 va- 
cancy in an initial bulk cell containing m chemical units 
of ABO3 relative to the energy of m — 1 chemical units 
of ABO3 and one chemical unit of AO. 



Figure [12] presents Ey computed as described above, 
for m — rf (corresponding to a series of supercells of 
size rj = 4, 6, . . . , 14 primitive unit cells on a side, with 
between 317 and 13 717 atoms each), plotted as a func- 
tion of l/f] 2 . Following a least-squares fit, the unam- 
biguous linear behavior in the plot indicates that the net 
interaction decays as l/i] 2 , with the ordinate intercept 
(I/77 2 — * 0, or 77 — > 00) providing the extracted value of 
Ey for the incorporation of excess AO as fully isolated 
BO2 vacancy complexes. The result within our model is 
Ey — 3.61 eV per excess AO unit, far exceeding the en- 
ergy of formation for even the least favorable RP phase 
configuration (162 meV per excess AO chemical unit for 
71 — * 00 in the [0 - ] configuration). The positive sign 
and large magnitude of this energy indicate that A-rich 
ABO3 will preferentially form planar AO stacking faults 
rather than isolated BO2 vacancy complexes. This result 
supports the present work's focus on such stacking faults 
and RP phases. 



VI. SUMMARY AND CONCLUSION 

This work presents a detailed and exhaustive study 
of the zero-temperature structures of the homologous 
A„+iB n 3 „ + i series of Ruddlesden-Popper (RP) phases 
in a model antiferrodistortive perovskite for the great- 
est range of n considered to date. We also consider, for 
a specific Glazer system, all possible octahedral orienta- 
tions within RP phases of the material. Our work in- 
troduces and follows a general program which may be 
replicated for RP phases of perovskites within any Glazer 
system and applied in both ab initio and empirical po- 
tential studies. 

We begin our particular considerations with a care- 
ful enumeration and symmetry analysis for all thirty-two 
combinatorially possible configurations of octahedral ro- 
tations within RP phases of an a~a~a~ Glazer system. 
Symmetry arguments reduce these thirty-two possibili- 
ties to five distinct configurations, which we denote as 
[0+], [0-], [§+], [7T+], and [tt-]. 

We then proceed to determine, within a shell-potential 
model, the ground-state structures of these five symme- 
try-distinct configurations for all values of n = 1 ... 30. 
We find that, for our specific model, only four config- 
urations are actually distinct, as one of the configura- 
tions ([f + ]) is unstable to immediate relaxation to an- 
other configuration {[it~ ]). After additional perturbation 
and quenching tests, we identify three of these remain- 
ing configurations as stable [0 + ], [7r + ]), and one as 
meta-stable ([0 - ]). In agreement with experimental ob- 
servation across a range of perovskites, we find that, for 
all values of n, all of the configurations of RP phases pro- 
vide a much more stable option for the incorporation of 
excess A-oxide than the formation of point defects, such 
as BO2 vacancy complexes. We also find that, within 
our model, the two most stable configurations ([ir~] and 
[0 + ]) lead to stacking faults which are energetically pre- 
ferred over phase separation into the bulk perovskite and 
bulk A-oxide. Both of these findings are consistent with 
the experimental observation of RP phases in a wide vari- 
ety of antiferrodistortive perovskites, including titanates, 
ruthcnates, manganites, and niobates. 

In terms of the interactions between stacking faults, we 
find that the interaction varies inversely with the distance 
between excess AO planes, a result which we believe to 
be generic, at least for the a~a~a~ system. For our 
particular model perovskite, the form of the interaction 
between stacking faults is attractive. In perovskites with 
such an attractive interaction, we find that ABO3 with 
relatively low stoichiometric excess of species A mani- 
fests intergrowths between the n = 2 RP phase and bulk 
perovskite. When such intergrowths are observed experi- 
mentally, the interaction is likely attractive in correspon- 
dence with our model material. 

Additionally, we observe that the five symmetry-dis- 
tinct configurations have significant differences in ener- 
gies of formation: ~30 meV per rotational degree of free- 
dom for our model perovskite, appreciably greater than 
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the energy scales associated with orientational disorder 
in our model perovskite, 4-16 meV. This leads to an in- 
triguing conjecture. For low densities of stacking faults 
within some perovskitcs, for a range of temperatures 
above their orientational phase transitions, the octahe- 
dra in bulk regions would exhibit randomized behavior 
about a high-symmetry mean configuration while octa- 
hedra near stacking faults would exhibit more ordered 
behavior about a lower-symmetry mean. Alternatively 
for higher densities of such faults, this increase in energy 
scale would correlate to an increase of the transition tem- 
peratures associated with octahedral rotations. These 
effects may be particularly noticeable in materials with 
higher orientational energy scales than our model per- 
ovskite, such as lanthanum aluminate with its transition 
temperature of 800 K. 

Finally, we propose a simple physical mechanism to 
explain the strong dependence of the interfacial energy 



on the rotational state of the octahedra: some configu- 
rations result in movement of neighboring like-charged 
oxygen ions directly toward each other and thus are high 
in energy, whereas others result in movement of oxygen 
ions past each other and thus are low in energy. 
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